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Mott transition and ferrimagnetism are studied in the Hubbard model on the anisotropic kagome 
lattice using the variational cluster approximation and the phase diagram at zero temperature and 
half-filling is determined. The ferrimagnetic phase rapidly grows as the geometric frustration is 
relaxed, and the Mott insulator phase disappears in moderately frustrated region, showing that the 
ferrimagnetic fluctuations stemming from the relaxation of the geometric frustration is enhanced by 
the electron correlations. In metallic phase, heavy fermion behavior is observed and mass enhance- 
ment factor is computed. Enhancement of effective spatial anisotropy by the electron correlations is 
also confirmed in moderately frustrated region, and its effect on heavy fermion behavior is examined. 
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Electron systems with geometric frustration play an 
important role to reveal the nature of the Mott trantision 
and physics related to it, e.g., heavy fermion behavior 
and spin liquids. When spatial anisotropy is introduced 
in systems with geometric frustration, the interplay be- 
tween the spin fluctuations and Mott transition appears 
as a new feature and provides unique phenomena which 
take place neither in the unfrustrarcd nor fully frus- 
trated systems. A reentrant behavior of the Mott transi- 
tion observed in the k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 un- 
der pressure^ is an interesting example realized on an 
anisotropic triangular-lattice, where that behavior stems 
from the enhancement of the antifcrromagnetic fluctua- 
tions due to the electron correlations^ 

As for the kagome lattice, which is a prototype of frus- 
trated systems, a fully frustrated case has been theo- 
retically studied in detail^— however the issues related 
to the anisotropy have been considered only recently. 
Spin correlations and magnetic properties in the metallic 
phase, as well as the Mott transition, have been ana- 
lyzed using the cellular dynamical mean field theory^, 
where the Mott transition point was determined, and 




FIG. 1: (Color online) (a) Anisotropic kagome lattice. 
In our lattice geometry, the three sites 1,2, and 3 form 
an equilateral triangle of the unit length and the dashed 
lines are along the x direction. Inside the hexagon (dotted 
line) and square (dash-dotted line) are the 12- and 6-site 
clusters, respectively, which will be used in our analysis, 
(b) The first Brillouin zone of the anisotropic kagome lattice. 



enhancement of spatial anisotropy and spin correlations 
were observed in the metallic phase. Such enhancement 
suggests that the spin correlations may be also enhanced 
above the Mott transition, giving rise to the extension of 
the ordered (ferrimagnetic) phase. Therefore, if the Mott 
transition itself persists without being veiled by the ferri- 
magnetic phase remains to be examined. Also, the effect 
of the enhanced anisotropy on the heavy fermion behav- 
ior is worth being studied. 

In this paper, we investigate the ferrimagnetism and 
Mott transition on the anisotropic kagome lattice using 
the variational cluster approximation (VCA)2r— , which 
is formulated based on a rigorous variational principle 
and exactly takes into account the short-range correla- 
tions. We determine the phase diagram at zero tem- 
perature and half-filling. We show that, in moderately 
frustrated region, the ferrimagnetic phase rapidly grows 
down to the metal-insulator phase boundary, indicating 
that the spin correlations stemming from the relaxation 
of the frustration is enhanced by the electron correlations 
and the Mott insulator (MI) phase disappears. In the 
metallic phase, heavy fermion behavior is observed and 
the mass enhancement of the quasiparticle is computed. 
Effective spatial anisotropy becomes also larger due to 
the electron correlations, in agreement with the previous 
study4 This effect gives rise to an enhancement of the 
anisotropy of the effective masses of the quasiparticles. 

The Hamiltonian of the Hubbard model on the 
anisotropic kagome lattice (see Fig. Q} reads 



H — tijcl a Cj a 



(i) 



where tij = t between the sites 1 and 2, 3 and tij = t! 
between the sites 2 and 3, U is the on-site Coulomb repul- 
sion, and fi is the chemical potential. The annihilation 
(creation) operator for an electron at site i with spin a is 

denoted as Cj„ (c| CT ) and = c\ a Cia- The system corre- 
sponds to the fully frustrated kagome lattice at t' /t = 1, 
and frustration becomes weaker with decreasing t' /t. The 
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end member att'/t = is a decorated square lattice. The 
energy unit is set as t = 1 hereafter. 

We use VCA^r— to examine the phase diagram and 
behavior of the quasiparticles in the metallic phase at 
zero temperature. VCA is an extension of the cluster 
perturbation theory^ based on the self-energy-functional 
approach*^ This approach uses the rigorous variational 
principle <5fit[E]/(5E = for the thermodynamic grand- 
potential fit written as a functional of the self-energy E 



O t [E] =F[E]+Trln(-(G 1 -E)- 1 ). 



(2) 



In the above expression, F[Y,] is the Legendre transform 
of the Luttinger-Ward functional^ and the index t de- 
notes the explicit dependence of fit on all the one-body 
operators in the Hamiltonian. The stationary condition 
for fit [E] leads to the Dyson's equation. All Hamiltonians 
with the same interaction part share the same functional 
form of F[S\, and using that property F[S] can be eval- 
uated from the exact solution of a simpler Hamiltonian 
H', though the space of the self-energies where F[H] is 
evaluated is now restricted to that of H' . In VCA, one 
uses for H 1 a Hamiltonian formed of clusters that are dis- 
connected by removing hopping terms between identical 
clusters that tile the infinite lattice. A possible symme- 
try breaking is investigated by including in H' the cor- 
responding Weiss field that will be determined by mini- 
mizing the grand-potential fit. Rewriting F[S\ in terms 
of the grand-potential fi' = fi^[E] and Green function 
G' _1 = Gq _1 — E of the cluster Hamiltonian H' , the 
grand-potential is expressed as 

fit(t') = n'-f pe 5 " V bidet (1 + (Go" 1 - G'^G') 
Jc 2tt K 

(3) 

and is now a function of t'. The functional trace has be- 
come an integral over the diagonal variables (frequency 
and superlattice wave vectors) of the logarithm of a deter- 
minant over intra-cluster indices. The frequency integral 
is carried along the imaginary axis and 5 — > +0. The 
stationary solution of fit(t') and the exact self-energy of 
H' at the stationary point, denoted as E*, are the ap- 
proximate grand-potential and self-energy of H in VCA, 
and physical quantities, such as expectation values of the 
one-body operators, are calculated using the Green func- 
tion Go -1 — E*. In VCA, the restriction of the space of 
the self-energies E into that of H' is the only approxi- 
mation involved and short-range correlations within the 
cluster is exactly taken into account by exactly solving 
H'. 

In our analysis, the 6- and 12-site clusters in Fig. [TJa) 
are used to set up the cluster Hamiltonian H' , and to 
study the ferrimagnetism, the Weiss field 



(4) 



with sign(i) = —1 for the site 1 and sign(i) = 1 for 
the sites 2 and 3, is also included. In the stationary 
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FIG. 2: (Color online) Phase diagram of the Hubbard model 
on the anisotropic kagome lattice at zero temperature and 
half-filling as a function of t' and U obtained by VCA, where 
the 12-site cluster is used. Lines are guides to the eye. 

point search of f2(/z', /if), which we denote as the grand- 
potential per site, the Weiss field /if and the cluster 
chemical potential \j! are treated as the variational pa- 
rameters, where the latter should be included for the 
thermodynamic consistency^ During the search, the 
chemical potential of the system [i is also adjusted so 
that the electron density n is equal to 1 within 0.1%. In 
general, a stationary solution with Zip 7^ corresponding 
to the ferrimagnctic state and that with h-p — corre- 
sponding to the paramagnetic state are obtained, and the 
ground-state energies per site fi + /m are compared to de- 
termine which solution (ferrimagnetic or paramagnetic) 
is stable. The density of state 
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is also calculated to examine the gap. To be precise, 
D(u>) is calculated for 77 = 0.2, 0.1, and 0.05, and rj — >• 
limit is evaluated by the standard extrapolation method. 
The numerical error after this extrapolation is estimated 
to be of order 10~ 3 , so the gap is identified as the region 
of us around ui ~ where the extrapolated D(uj) is less 
than 10~ 2 . 

In Fig. [2 we show the phase diagram at zero temper- 
ature and half-filling obtained from this analysis using 
the 12-site cluster. The critical interaction strength Up 
separating the ferrimagnetic and MI phases rapidly de- 
creases in the moderately frustrated region t' = 0.5 ~ 0.7, 
showing that the ferrimagnetic fluctuations due to the re- 
laxation of the geometric frustration is enhanced by the 
electron correlations. In this region of t' , the ferrimag- 
netic phase is an insulator since there is a gap, and the 
transition between the ferrimagnetic and paramagnetic 
(including MI) phases is a level crossing because the fer- 
rimagnetic solutions exist also U < Uf even though it 
is energetically disfavored there. The critical interaction 
strength Umi separating the MI and metallic phases is 
slightly smaller than the noninteracting band width W, 
where W = 6 at t' = 1 and W = 4v^ ~ 5.66 at t' = 0. 
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FIG. 3: (Color online) Mott gap A as a function of U at (a) 
t' = 1.0, (b) t' = 0.8, and (c), (d) t' = 0.6. For t' = 0.6, order 
parameter M is also included. The vertical lines separate the 
ferrimagnetic and MI phase. 

Umi decreases as the geometric frustration is relaxed, and 
the slope becomes steeper in moderately frustrated re- 
gion. At t' = 0.5, U F = 4.0 while U M i = 4.1 so the MI 
phase has disappeared. Taking into account the drastic 
growth of ferrimagnetic phase and the fact that W re- 
mains almost the same, the decrease of Umi according to 
the relaxation stems from the ferrimagnetic fluctuations. 

In Fig. [3] we show the Mott gap A and ferrimagnetic 
order parameter M as functions of U for t' = 1.0, 0.8, 

and 0.6, where M = X)a=i(< "at > — < n ai >) an d 
< n aa > is the expectation value of n a(T with a =1, 2, and 
3 being the sites in Fig. HJa). The results obtained using 
the 6-site cluster are also given to see the cluster size de- 
pendence quantitatively. A monotonically decreases as 
U decreases. Umi is always larger in the 12-site analy- 
sis, but the differences are less than 20% of W. The gap 
Af in the ferrimagnetic phase is slightly larger than A 
if U is the same. For example, for U = 10, A = 4.35 at 
t' = 0.8 while A F = 4.91 at t' = 0.7 and A F = 5.54 at 
t' = 0.5 in the 12-site analysis. Up is almost the same for 
the 12- and 6-site clusters at t' = 0.6 and M is always 
smaller for the 12-site cluster. In general, Umi is larger 
for larger clusters, since the kinetic energy of the clus- 
ter Hamiltonian can be larger for larger clusters. As for 
Up, when spin correlations are highly suppressed due to 
the frustration, cluster wave functions with small ferri- 
magnetic fluctuations play an important role to examine 
near the true minimum of the effective potential, so Up 
is smaller for larger clusters. When the geometric frus- 
tration is moderate and spin correlations are not largely 
suppressed, the difference of cluster kinetic energies due 
to the cluster size becomes more important to determine 
the phase boundary, so Up becomes larger for larger clus- 
ters. U F = 9.5 at t' = 0.7 and U F = 3.6 at t' = 0.5 for 
the 6-site cluster, which is consistent with this general ar- 
gument on the cluster size dependence. Relatively large 
difference of U-p in strongly frustrated region is due to 
strong suppression of the spin correlations. 

In Fig. |||we show the spectral weight function p(u>, k) 
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FIG. 4: Spectral density at t' = 0.75 for (a), (b) U = 8 
(ferrimagnetic state), (c) U = 6 (MI state), and (d) U — 3 
(metallic state) along the dotted line in Fig. \Hb). The 
Lorentzian broadening with r\ — 0.15t is used in all the 
cases. In (a) and (b), the solid lines are the mean-field 
SDW dispersion for the same values of U, t' and /j,. In 
(d), the solid lines are the noninteracting band structure. 
In (a), (b), and (c) the peaks are scaled by 5 compared to (d). 



calculated using the 12-site cluster for solutions corre- 
sponding to (a), (b) the ferrimagnetic phase (up and 
down spin parts are plotted separately), (c) the MI 
phase, and (d) the metallic phase at t' = 0.75. In (a) 
and (b), the mean-field spin-density- wave (SDW) disper- 
sion is also included (solid lines) to see its general fea- 
tures, where M = 0.92 in the mean- field solution while 
AI = 0.72 in VCA. In (d), the noninteracting band struc- 
ture is also plotted (solid lines) for comparison. In (a) 
and (b), the dispersion is largely affected due to the elec- 
tron correlations compared to the mean-field solution. In 

(c) , the SDW dispersion disappears and the spectral func- 
tion displays a Mott gap across all wave vectors. The gap 
is smaller compared to (a) and (b) since U is smaller. In 

(d) , we notice that the interacting bands slightly shrink 
toward the Fermi surface, leading to heavy fermion be- 
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FIG. 5: (Color online) (a) Ratio r as a function of £7 for 
t' = 1.0, 0.8, and 0.6. The two arrows indicate the value 
of r for noninteracting band at t' — 0.5 (r = 1.24) and at 
t' = 0.4 (r = 1.30). (b)~(d) Mass enhancement factor m*/m 
in the x and y directions as functions of U for (b) t' = 1.0, 
(c) t' = 0.8, and (d) t' = 0.6. The 12-site cluster is used. 

havior. 

To study it in detail, we compute the mass enhance- 
ment factor m*/m along the x and y directions in k space, 
where the x direction corresponds to the direction of t' 
hopping in real space. Near the Fermi surface the posi- 
tion of the peak (w, k) = (u>p + Su>, kp + 5k) of p(ui, k) 
changes according to the relation \8uj\ = kp/m*8k and 
the band mass m* is calculated using this relation. We 
also compute the ratio of the Fermi momenta in the x and 
y directions, r = k y p/k x p. As t' decreases, the nonin- 
teracting Fermi surface slightly shrinks in the x direction 
and slightly evolves in the y direction^, so r is a measure 
of the anisotropy including the effect of the electron cor- 
relations. Even though precise values of these quantities 
may depend on a lattice geometry, we show in Fig. [5] r 
and m* /m in the x and y directions as functions of U 
for t' = 1.0, 0.8, and 0.6 in our lattice geometry, to see 
general features about the effect of the electron correla- 
tions on these quantities. When U ~ 3, r rapidly grows 



for t' = 0.6, indicating that the effective anisotropy is 
enhanced due to the electron correlations in moderately 
frustrated region. In fact, the value of r at t' = 0.6 
and U = 3.5 is equal to that of noninteracting band at 
t 1 = 0.4. This enhancement of effective anisotropy, re- 
ported also in Ref. 8, indicates that the spin fluctuations 
due to the relaxation of the geometric frustration are also 
enhanced by the electron correlations. Within our anal- 
ysis, this indication is consistent with the rapid growth 
of the ferrimagnctic phase above Umi, and it is demon- 
strated by the analysis of the spin correlations in the 
metallic region^ As is shown in Fig.[SJb)~(d), the heavy 
fermion behaviors are observed in all cases. The growth 
of r affects also m* /m since kp enters into the calcula- 
tion of in* . In general, in* is enhanced in the y direction 
and suppressed in the x direction since the Fermi surface 
shrinks in the x direction and evolves in the y direction. 
This appears largely in moderately frustrated region due 
to the rapid growth of r, as is observed in for U ~ 3.5 
in Fig. [5jd). Therefore the anisotropy of the effective 
masses is enhanced in moderately frustrated region. 

In summary we have investigated the ferrimagnetism 
and Mott transition on the anisotropic the kagome lat- 
tice using VCA. The phase diagram at zero temperature 
and half-filling is determined. The ferrimagnetic phase 
rapidly grows in moderately frustrated region and the 
MI phase disappears there. In the metallic phase, heavy 
fermion behavior is studied and the mass enhancement is 
computed. Enhancement of spatial anisotropy due to the 
electron correlations is also observed for moderately frus- 
trated region and its effect on the heavy fermion behavior 
is discussed. Thus, the interplay between the spin cor- 
relations and Mott transition are quantitatively studied 
above and below the metal-insulator transition. 
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